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摘要 明 安 图 射电 频谱 日 像 仪 (Mingantu Spectral Radioheliograph, MUSER) 能 够 在 0.4-15 GHz 超 宽频 带 内 

实现 高 时 间 、 高 空间 、 高 频率 分 辨 率 的 太阳 射电 成 像 . 而 射电 亮 温度 是 描述 太阳 物理 过 程 的 一 个 重要 的 参数 ， 

在 研究 不 同 射电 辐射 机 制 、 太 阳 磁 场 以 及 太阳 爆发 过 程 中 非 热 粒子 加 速 等 问题 上 有 着 非常 重要 的 作用 , 因此 必 

须 对 MUSER 观 测 的 图 像 进 行 亮 温度 定 标 . 将 介绍 一 种 适用 于 射电 日 像 仪 图 像 强度 定 标的 方法 . 太阳 射电 图 像 中 

包含 着 太阳 圆 盘 的 结构 信息 , 利用 射电 日 像 仪 短 基线 的 可 视 度 函 数 拟 合 第 一 类 贝 塞 尔 函 数 , 可 以 得 到 图 像 中 宁 

静 太 阳 圆 盘 的 射电 半径 和 强度 , 再 利用 瑞 利 - 金 斯 定律 和 每 天 的 太阳 射电 流量 可 以 计算 得 到 每 天 图 像 的 定 标 因 

子 Ge, 从 而 实现 对 MUSER 图 像 强度 的 定 标 . 将 该 方法 应 用 到 MUSER 的 实际 观测 数据 中 , 包括 宁静 太阳 和 有 太 

阳 射 电 爆 发 等 不 同 的 情况 , G6。 的 误差 基本 不 超过 10%, 得 到 的 宁静 太阳 亮 温 度 与 其 他 宁静 太阳 的 结果 有 具有 和 较 高 

的 相关 性 , 表明 了 此 方法 的 可 行 性 和 有 效 性 . 

关键 词 太阳 : 射电 辐射 , 仪器 : 射电 日 像 仪 , 仪器 : 定 标 

中 图 分 类 号 : P162; ”文献 标识 码 : A 
1 5| 言 tion, CME) 以 及 各 种 热 和 非 热 过 程 . 太阳 爆发 被 认 

太阳 射电 观测 可 以 提供 从 太阳 色 球 到 日 地 空 为 是 由 于 太阳 磁场 的 拓扑 重组 或 磁场 重 联 导 致 的 
间 广 阔 区 域 中 有 关 等 离子 体 和 高 能 粒子 动力 学 行 瞬时 能 量 释 放 过 程 回 . X45] £X XL Ju ze HH, 非 热 粒 子 
为 等 信息 , 这 是 其 他 观测 手段 所 不 具备 的 , Dj, —— C EEROROION EREK, EA PEDERI SURG BT 
射电 观测 是 研究 太阳 物理 的 重要 探测 手段 站. 太阳 。” 段 存 在 大 量 的 耀 斑 能 量 积聚 因此 , 对 不 同类 型 
射电 观测 范围 包括 宁静 太阳 、 耀 班 、 日 冕 物质 抛 ”太阳 射电 暴 的 观测 对 于 理解 粒子 加 速 的 机 制 以 及 
射 和 行星 际 粒子 等 , 而 太阳 射电 辐射 由 热 粒子 ( 自 ”太阳 爆发 的 起 源 和 演化 过 程 至 关 重 要 四 . 太阳 射电 
由 -自由 辐射 )、 中 等 相对 从 粒子 (回旋 同步 加 速 辐 观测 至 今 已 有 七 十 多 年 的 历史 , 观测 设备 由 早期 
射 ) 以 及 不 稳定 的 各 向 异性 粒子 (如 电子 束 产 生 等 离 。 ”的 射电 流量 计 、 冉 电 频 谱 仪 到 现在 的 射电 频谱 日 
子 体 辐射 、 损 失 锥 产生 电子 回旋 脉 泽 辐射 ) 产 生 癌 . 像 仪 . 明 安 图 射电 频谱 日 像 仪 (Mingantu Spectral 
此 外 , 太阳 射电 爆发 包含 的 信息 与 爆发 类 型 有 关 ， Radioheliograph, MUSER) 是 太阳 专用 的 射电 望 远 
WARNER H2% jhi (Coronal Mass Ejec- 镜 , 采用 综合 孔径 成 像 技术 , 在 0.4-15 GHz 超 宽频 
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65 卷 X x 学 d 
带 内 同时 实现 了 时 间 、 空 间 、 频 率 的 高 分 辩 率 太 位 校准 
EE X 8 671, 宽带 频谱 射电 成 像 可 以 提供 关于 太阳 中 的 两 个 20 m 的 低频 天 线 和 一 个 16 mit 4 i 
日 冕 磁场 、 能 量 释放 机 制 、 高 能 粒子 的 产生 和 传 收 机 的 高 频 天 线 ; 
播 以 及 相应 的 辐射 机 制 等 方面 的 丰富 信息 , 并 可 于 MUSER 相 关 接 
以 进一步 研究 和 揭示 太阳 色 球 层 和 日 冕 的 辐射 机 关 , 系统 增益 对 相关 结果 的 影响 不 大 , 因此 本 研 
EE. 而 射电 亮 温度 是 研究 太阳 物理 过 程 的 一 个 非 中 并 没有 考虑 系统 增益 的 校准 . 
常 重要 的 物理 参数 . 因此 , 对 射电 日 像 仪 的 图 像 进 当前 MUSER 已 经 
行 亮 温度 定 标 是 非常 必要 的 . 爆发 事件 , 包括 太 
MUSER 采 用 综合 孔径 成 像 方式 对 太阳 在 射电 射电 爆发 事件 


波段 上 成 像 ， 当 观 测 视 场 很 小 时 , 其 观测 天 区 的 
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[310 此 外 , MUSER 也 可 以 利用 正在 建设 


行 MUSER 阵 列 的 校准 


| 冷 接 


7,11] 由 


改 机 对 信号 进行 2 比特 量化 和 相 
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观测 到 很 多 类 型 的 太 B 
阳 泡 斑 之 前 和 期 间 的 准 周期 脉动 
5 、Fibers 爆 发 结构 、Spike 爆 发 结构 
以 及 Lace 爆 发 结构 等 03. 然而 , 当 


ET 


a 


前 MUSER 的 图 像 


射电 图 像 T(,) 可 通过 以 下 的 二 维 传 里 叶 变 换 得 还 没有 进行 强度 定 标 , 为 了 充分 利用 MUSER 图 像 
$i: 在 太阳 物理 研究 中 的 价值 , 因此 , 本 文 针 对 MUSER, 
| 的 图 像 提 出 了 一 种 强度 定 标 方法 . 
Pam = ff Ve udv, 在 太阳 射电 图 像 的 亮 温度 定 标 研究 中 , 西伯 利 
(1) 亚 射 电 日 像 仪 (Siberian Radioheliograph, SRH) 通 
过 分 析 太 阳 脏 图 的 强度 直方 图 , 获取 直方 图 中 包 
其 中 V(w,w) 为 可 视 度 函 数 ， u 和 wv 分 别 是 垂直 于 相位 4 含 的 两 个 峰值 ， 它们 分 Y IDEA 应 着 Mim 5 m3 cris 
参考 位 置 方向 的 平面 的 坐标 轴 分 量 , v 是 朝 着 通过 DE TOR 1 
EN N : 太阳 圆 盘 的 强度 值 (任意 单位 ). 再 利用 Zirin 等 人 
原点 的 平面 定义 的 北方 , v 是 朝向 东方 , I 和 m 是 对 在 1986_1987 年 太阳 黑子 极 小 期 通过 单口 径 天 线 
m Mich E E : E 年 黑子 极 小 期 通过 FRZ 
应 着 图 像 平面 中 的 两 个 坐标 轴 分 量 , uvh "PPS " Wl DESENES 
D 、 i 、 观测 得 到 的 宁静 太阳 的 亮 温 度 值 9, 通过 以 下 公式 
HRI, 7 为 虚数 单位 ,可 视 度 函 数 Yu 通过 阵列 ree nag rt 
中 任意 两 天 线 的 信号 相关 得 到 , 其 对 应 观测 天 区 射 
电 图 像 T(1,m) 的 侍 里 叶 成 分 ; 把 射电 阵列 的 所 有 天 o= [S00 - 9e | op s), (3) 
线 两 两 相关 (天 线 对 ) 得 到 的 可 视 度 函 数 进行 傅 里 叶 Qiss(v) - Qv) 
反 变 换 即 可 得 到 观测 天 区 的 射电 图 像 . 由 于 实际 中 CC(w) 为 定 标 后 的 图 像 值 , > 表示 频率 , Tuas (v) 
中 射电 阵列 的 天 线 个 数 有 限 , 不 能 把 射电 图 像 的 。 ”是 字 静 太阳 的 亮 温度 是 通过 Zirin 等 9 观测 到 的 
所 有 傅 里 叶 成 分 都 得 到 . 在 这 种 情况 下 , 缺失 了 部 结果 得 到 , 并 且 每 个 频率 只 有 一 个 固定 值 , Qus (v) 
分 侍 里 叶 成 分 的 射电 图 像 TD(1,m) 被 称 为 脏 图. 因 和 Qieva(w) 分 别 为 强度 直方 图 中 的 两 个 峰值 . 
此 , 需要 进一步 做 去 卷 积 处 理 得 到 真实 的 射电 图 针对 射电 源 的 强度 定 标 , 扩展 欧文 斯 谷 太阳 阵 
P e ierat : 列 (Expanded Owens Solar Array, EOVSA) 
然而 ， 在 实际 观测 中 ， (1) 式 中 的 一 对 基线 (天 线 可 以 通过 射电 源 的 流量 E 密度 和 射电 源 的 面积 得 到 
p 和 天 线 g) 的 可 视 度 函数 可 表示 如 下 : 对 应 射电 源 的 亮 温 度 [19. 


pq 一 


r? = | VPI] g”? exp iloi + p30)], (2) 


然而 ，SRH (Statistical Result Histogram) 的 


定 标 方法 只 适合 宁静 太阳 和 微弱 爆发 事件 的 图 像 
其 中 rz 表示 观测 得 到 的 基线 (z 一 g) 的 干涉 响应 ， 定 标 , 并 且 太 阳 峰 年 和 谷 年 的 宁静 太阳 亮 温 度 是 不 
Vif n] 视 度 函数 的 幅度 , 92 是 干涉 元 的 幅度 增 一 样 的 . 因此 , 基于 每 天 的 太阳 射电 流量 , 本 文 提出 
益 , yp 名 是 可 视 度 函数 的 相位 ,p34 是 干涉 元 仪器 一 种 适用 于 射电 日 像 仪 图 像 亮 温度 定 标 的 方法 , 在 
效应 带 来 的 相位 误差 . 因此 , 在 成 像 之 前 需要 对 第 2 节 将 详细 介绍 定 标的 方法 和 原理 , 第 3 节 对 实际 
可 视 度 函 数 进行 校准 . 当前 MUSER 利 用 风云 系列 的 MUSER 观 测 数据 进行 定 标 , 展示 定 标 后 的 结果 


同步 卫星 以 及 太阳 圆 盘 上 射电 爆发 源 来 进行 相 
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并 加 以 分 析 , 最 后 为 结论 . 


65 卷 周志 超 等 : 明 安 图 射电 频 


A 
2.1 


图 像 强度 定 标 方 法 

宁静 太阳 圆 盘 的 射电 半径 和 强度 

基于 EOVSA 预 估 射 电 图 像 中 宁静 太阳 圆 盘 的 
射电 半径 和 强度 的 思路 !, 射电 日 像 仪 观测 的 太阳 
图 像 数 据 中 包含 了 太阳 圆 盘 的 结构 信息 , 阵列 中 短 
基线 的 可 视 度 函数 反映 的 是 图 像 的 大 尺度 结构 . 
此 , 短 基线 的 可 视 度 函数 应 该 符合 圆 盘 的 傅 里 叶 变 
换 规律 , 满足 第 一 类 贝 塞 尔 函 数 , 即 幅 度 具 有 以 下 
关系 : 


2 ita) | (4) 


a=Q (532 

z 
KENE —3 BURG BUG Qo MA 
值 , > 是 和 太阳 圆 盘 的 直径 有 关 的 量 , 表达 式 为 : 


(5) 


其 中 的 轧 为 弧度 制 下 的 太阳 直径 , p 是 波长 和 单位 下 
的 wv 距离 , Bo = vu 十 v2, DD 不 是 太阳 的 光 球 层 
直径 , 而 是 不 同 频率 的 射电 直径 . 通过 调整 (4) 式 和 
(5) 式 中 圆 盘 的 强度 Qo 和 直径 万 可 以 调整 对 应 的 曲 
线 4, 利用 实际 的 短 基线 可 视 度数 据 拟 合 出 一 条 最 
接近 的 曲线 4, 从 而 得 到 对 应 的 圆 盘 强度 和 直径 . 
SRH 通 过 分 析 射 电 日 像 仪 脏 图 的 直方 图 , 分 
别 统计 图 像 中 1.2 倍 太阳 半径 以 外 和 0.8 倍 太阳 半 
径 以 内 区 域 的 强度 值 , 这 两 部 分 的 峰值 分 别 对 应 
着 图 像 中 天 空 背景 强度 值 和 宁静 太阳 圆 盘 强 度 值 
TuuistGei,， 其 中 Ge 为 射电 日 像 仪 图 像 的 定 标 因 子 . 
但 实际 由 于 宁静 太阳 圆 盘 的 强度 是 天 空 背景 的 上 
万 倍 , 天 空 背景 的 强度 不 能 很 灵敏 地 体现 在 射电 日 
像 仪 图 像 中 . 因此 , 在 利用 SRH 的 方法 时 我 们 只 对 
MUSER 图 像 中 0.8 倍 太阳 半径 以 内 的 区 域 进行 了 
统计 , 如 图 1 所 示 . 图 1 (a) 为 2016 年 7 月 5 日 MUSER 
在 1.7125 GHz 未 定 标 的 脏 图 I?(1,m), 其 中 黄 线 以 
内 表示 光学 太阳 圆 盘 区 域 , 红色 虚线 则 表示 0.8 倍 
太阳 半径 , 颜色 棒 表 示 其 脏 图 中 数值 大 小 (任意 单 
位 ), 图 1 (b) 为 对 图 1 (a) 中 0.8 倍 太阳 半径 以 内 区 域 
的 强度 值 统计 结果 , 黑色 曲线 表示 统计 结果 , 红色 
虚线 为 拟 合 的 高 斯 曲线 ,高 斯 拟 合 的 期 望 值 对 应 


z—mDp, 


3 HH 


像 仪 图 像 的 强度 定 标 方法 


着 第 2 个 峰值 Qievz, 即 通 过 统计 直方 图 的 方法 得 到 
Tue Gc = 774806 + 152638. 


a (b): 
o d 
E E 
E 
总 E 
p E 
2 E 
E 

J 当 
E j 

1.2x108 
Q (arbitrary unit) 

图 1 MUSER 在 1.7125 GHz 未 定 标 的 脏 图 I?(1, m) 和 其 强度 Q@ 的 直 


方 图 统计 结果 , 都 为 任意 单位 . 


Fig.1 Uncalibrated dirty image IP (l, m) of MUSER at 
1.7125 GHz, and the histogram statistical result of its 


intensity Q, both are in arbitrary unit. 


此 外 , 利用 拟 合 第 一 类 贝 塞 尔 函 数 的 方法 ， 
我 们 对 MUSER 在 2016 年 7 月 5 日 几 个 小 时 短 基线 观 
测 到 的 可 视 度 函数 幅度 进行 分 析 . 通过 改变 (4) 式 
和 (5) 式 中 圆 盘 的 强度 Co 和 直径 万 , 拟 合 出 一 条 最 
接近 观测 数据 的 曲线 4, 即 找 出 一 组 (Qo,D), 使 
得 其 差 的 均 方 根 最 小 , iR (Qo, D) 对 应 MUSER 图 像 
中 宁静 太阳 圆 盘 强度 值 TwiwGe 和 太阳 圆 盘 直 径 
D 图 2 为 MUSER 在 2016 年 7 月 5 日 几 个 小 时 短 基 
线 观 测 到 的 可 视 度 函 数 幅 度 以 及 拟 合 的 情况 , 其 中 
黑色 星 号 表示 MUSER 在 1.7125 GHz 的 观测 数据 ， 


lhttp://www.ovsa.njit.edu/wiki/index.php/Full_Disk_Simulations 
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65 卷 


红色 的 曲线 则 是 拟 合 得 到 的 最 佳 曲线 4, 其 对 应 圆 
iit ELTE D7335.2', TiuietGe = 199540 + 8961.2. 


期 


镜 网 络 (Radio Solar Telescope Network, RSTN)? 
每 天 的 数据 得 到 . 


对 比 上 面 两 种 求 宁静 太阳 圆 盘 强 度 值 TuuetG。 但 是 太阳 射电 图 像 并 不 是 一 直 都 是 一 个 平坦 
的 结果 , 可 以 发 现 , 通过 统计 图 像 强度 的 直方 图 得 的 面 源 , 圆 面 上 的 活动 区 会 导致 误差 , 太阳 射电 流 
到 的 宁静 太阳 圆 盘 强 度 值 偏 大 . 这 主要 原因 是 , 尽 量 主要 由 平坦 的 圆 面 和 活动 区 两 部 分 页 献 . 因此 ， 
管 2016 年 7 月 5 日 那天 太阳 相对 宁静 , 没有 射电 爆发 若 减 去 活动 区 页 献 的 射电 流量 , 则 利用 瑞 利 - 金 斯 
事件 , 但 是 太阳 圆 面 上 有 若干 个 活动 区 , 同时 , 从 射 公式 可 近似 得 到 太阳 圆 面 的 射电 亮 温 度 . 这 里 通 


电 图 像 上 可 以 发 现 太阳 圆 面 上 有 一 些 较 强 的 结构 
信息 , 这 样 就 导致 了 通过 统计 得 到 的 太阳 圆 盘 强 度 


过 射电 图 像 中 活动 区 的 强度 和 面积 以 及 整个 圆 面 
的 强度 , 可 以 计算 出 活动 区 贡献 的 射电 流量 值 . 由 


比 真 实 的 值 偏 大 . 相反 , 利用 短 基 线 可 视 度 函数 幅 


于 天 空 背景 辐射 的 亮 温度 约 为 3 KL, 远 远 小 于 太 


度 拟 合 第 一 类 贝 塞 尔 函数 曲线 的 方法 则 不 受 这 些 
因素 影响 , 拟 合 结果 只 由 图 像 中 包含 的 圆 盘 结构 信 
EE. 因此 , 这 种 方法 得 到 的 宁静 太阳 圆 盘 强度 
比 统计 直方 图 的 结果 更 准确 . 这 也 印证 了 SRH 通 过 
统计 图 像 强度 直方 图 的 定 标 方法 只 适用 于 宁静 太 
阳 以 及 相对 较 弱 射电 爆发 事件 图 像 的 定 标 [9. 


Visibilities Amplitude for Short Baselines 
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图 2 2016 年 7 月 5 日 MUSER 在 1.7125 GHz 的 可 视 度 函数 幅度 以 及 


拟 合 的 曲线 


Fig.2 Visibilities amplitude of MUSER at 1.7125 GHz on 
July 5, 2016 and curve obtained by fitting 


2.2 ”利用 太阳 射电 流量 来 进行 图 像 定 标 
RR 据 瑞 利 - 金 斯 定律 09, 可 以 得 到 太 B 
射电 流量 具有 以 下 关系 : 

e sS 
7 hn: (6) 
PTUACK BURGRUR, AHE, Ks rU HERE 
量 , 9 为 太阳 射电 半径 对 应 的 立体 角 ,5 为 太阳 的 射 
电流 量 . 而 太阳 射电 流量 5 可 以 通过 射电 太阳 望 远 


亮 温度 


EEN 


Ty 


?https:/ /www.ngdc.noaa.gov/stp/space-weather/solar-data/solar-features/solar-radio /rstn-1-second/ 
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阳 的 亮 温 , 为 了 简化 计算 不 考虑 天 空 背景 的 辐射 . 


TL HO 
MUSER 图 像 的 定 标 因子 Ge 可 以 表示 如 下 : 


Qu = T Ge E (7) 


区 域 


像 中 圆 面 上 平坦 


中 Qu 为 MUSER 未 定 标 图 


pi 


的 平均 强度 值 , 图 像 定 标 则 可 以 通过 除 以 定 标 因子 


3 ” 图像 的 定 标 结果 对 比 和 分 析 
3.1” 定 标 因 子 G 的 结果 分 析 

本 方法 中 的 定 标 因子 G. 是 未 定 标 图 像 的 强度 
与 射电 亮 温度 的 比值 , 在 理想 状态 下 , 定 标 因 子 应 
该 是 一 个 固定 的 值 . 然而 , 由 综合 孔径 的 原理 可 知 ， 
射电 阵列 只 是 对 射电 源 的 部 分 传 里 叶 成 分 进行 采 
FÉ, 且 对 于 同样 的 天 线 阵 列 , 采样 函数 会 随 着 太阳 
立 置 的 变化 而 变化 . 对 于 太阳 角度 较 高 的 时 间 段 
内 , 其 短 基线 的 uv 距离 变化 较 小 , 从 而 对 于 射电 图 
像 的 总 强度 随 着 时 间 变 化 的 幅度 较 小 , 因此 , 在 这 
一 段 时 间 内 的 定 标 因子 相对 比较 稳定 . 然而 , 在 实 
际 观测 过 程 中 , 由 于 天 线 伺服 故障 以 及 部 分 基线 的 
相关 结果 不 理想 , 导致 每 天 的 采样 函数 不 一 样 , 因 
此 , 需要 每 天 对 MUSER 的 图 像 进行 强度 定 标 . 

对 MUSER 在 2016 年 7 月 5 日 的 观测 数据 进行 求 
G。， 用 2.1 节 方法 拟 合 得 到 1.7125 GHz 的 太阳 射 
电 直 径 为 35.2, 根据 RSTN 数 据 对 应 1.7125 GHz 
的 太阳 射电 流量 为 57.77 sfu. 利用 2.2 节 中 的 定 标 方 
iE, 得 到 对 应 的 定 标 因子 G. 如 图 3 所 示 , 取 平 均值 得 
到 当天 的 定 标 因 子 为 C。 = 6.526 + 0.658. 图 4 展示 


65 卷 同志 超 等 : 明 安 图 射电 频谱 日 像 仪 图 像 的 强度 定 标 方法 3 期 
了 2015 年 11 月 22 日 不 同 频率 求 到 的 Ge 值 情况 , 对 应 ” 在 那天 基本 处 于 宁静 状态 , 天 线 A2 接 收 到 的 功率 值 
1.26 GHz, 1.36 GHz, 1.46 GHz, 1.66 GHz. 1.71 GHz (任意 单位 ) 和 定 标 后 的 太阳 圆 盘 平均 亮 温度 也 基本 
和 1.86 GHz 的 定 标 因子 Ge 分 别 为 7.066 + 0.407、 维持 不 变 . 
9.122 土 0.607、8.247 土 0.718、9.387 土 0.737、8.919 土 el OE (3) 
0.703. 6.577 + 0.568. 从 以 上 计算 的 G. 值 来 看 , 可 S | gxi 
以 发 现 其 G。 值 在 一 天 内 都 保持 相对 稳定 , GRAE E pfoeeeeeeeettetesestooosssesooosssoooees 
1096 EA P4, 这 也 间接 反映 本 定 标 方法 的 稳定 性 和 可 2 
靠 性 $ 5.0x10'r j 
ut ot 
14F g 2.0x 2 (bj 
12r E ka 1.5x102 上 J 
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03:18:24 04:18:24 05:18:24 06:18:24 07:18:24 图 5 2016 年 7 月 5 日 MUSPR 图 像 强度 定 标 后 在 1.7125 GHz 的 太阳 


Time (UT) 


图 3 2016 年 7 月 5 日 MUSER 在 1.7125 GHz 对 应 的 G。 


Fig.3 The G. of MUSER at 1.7125 GHz on July 5, 2016 
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图 4 2015 年 11 月 22 日 不 同 频率 MUSER 图 像 的 定 标 因子 G。 值 


Fig.4 Calibration factor Gc values of MUSER images at 


different frequencies on November 22, 2015 


3.2 “” 定 标 后 MUSER. 的 图 像 分 析 

利用 本 方法 分 别 对 MUSER 在 2016 年 7 月 5 日 和 
2015 年 11 月 22 日 的 观测 数据 进行 图 像 定 标 . 这 两 天 
分 别 对 应 着 宁静 太阳 和 有 太阳 射电 爆发 的 两 种 情 
况 . 其 中 2016 年 7 月 5 日 定 标 后 太阳 圆 盘 平均 亮 温 度 
与 天 线 A2 接 收 到 的 功率 值 对 比如 图 5 所 示 . 图 5 (a) 
为 2016 年 7 月 5 日 天 线 A2 接 收 到 的 功率 值 (任意 单 
位 ), 图 5 (bp) 为 定 标 后 的 太阳 圆 盘 平均 亮 温 度 . 太阳 


PE 


hos 


45 
HA FRE CK 2E RE DREZ E: (a) 为 天 线 A2 接 
改 到 的 功率 值 (任意 单位 )，(b) 为 图 像 定 标 后 的 太阳 圆 盘 平 均 亮 温度 . 


Fig.5 Comparison between the average brightness 
temperature of solar disk at 1.7125 GHz after MUSER image 
intensity calibration and the power value received by a single 
antenna on July 5, 2016: (a) power value (in arbitrary unit) 
received by antenna A2, (b) average brightness temperature 


of solar disk after image intensity calibration. 


对 于 2015 年 11 月 22 日 的 MUSER 数 据 , 当天 太 
阳 有 射电 爆发 事件 , 活动 区 AR12454 在 01:37-02:27 
UT 发 生 了 GOES (Geostationary Operational En- 
vironmental Satellites) $& X4] ££ C5.12K I] E PE, 在 
03:17-03:25 UT Æ T GOESKX J 2X C1.22 B] E 
PE, 05:31-05:38 UT 发 生 了 GOES 软 X 射 线 C5.6 
ZR I XE BE. 而 MUSER 在 04:51-04:53 UT 和 05:35- 
05:38 UT 期 间 观 测 到 太阳 射电 爆发 , 在 相应 时 
间 段 , RSTN 在 1415 MHz、2695 MHz、4995 MHz、 
8800MHz、15400 MHz 等 频 点 上 也 观测 到 射电 爆 
R. 定 标 后 活动 区 AR12454 对 应 的 射电 源 在 04:51- 
04:53 UT 爆发 过 程 中 的 亮 温 度 以 及 单 天 线 的 功率 
直 变 化 情况 如 图 6 所 示 . 图 6 (a) 为 2015 年 11 月 22 日 
天 线 A2 的 功率 值 (任意 单位 ), 图 6 (b) 为 定 标 后 图 像 
中 最 大 的 亮 温 度 值 在 射电 爆发 期 间 的 变化 情况 , 即 
活动 区 AR12454 对 应 的 射电 源 亮 温 度 变 化 情况 . 同 
时 , 也 计算 了 图 6 (a) 中 的 功率 数据 与 图 6 (b) 中 的 亮 


— 
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温度 数据 的 相关 系数 R, 结果 如 图 7 所 示 . 从 图 中 可 
以 明显 看 出 , 射电 源 在 射电 爆发 期 间 的 亮 温 度 和 单 
天 线 接收 到 的 功率 值 (任意 单位 ) 具 有 相同 的 变化 趋 
A, 且 它 们 之 间 的 相关 系数 大 于 0.95. 
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图 6 2015 年 11 月 22 日 MUSER 图 像 强度 定 标 后 不 同 频率 射电 源 在 爆 


发 期 间 的 亮 温度 变 
(b) 图 


Fig.6 Brightness temperature change of radio source during 


化 : (a) 为 天 线 A2 接 收 到 的 功率 值 (任意 单位 )， 
像 定 标 后 射电 源 亮 温度 的 变化 . 


the burst period after MUSER image intensity calibration at 
different frequencies on November 22, 2015: (a) power value 
(in arbitrary unit) received by antenna A2, (b) brightness 
temperature change of radio source after image intensity 


calibration. 
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7 2015 年 11 月 22 日 不 同 频率 接收 到 的 功率 值 ( 任 ; 
度 定 标 后 射电 源 的 亮 温度 的 相关 性 


D 


图 


gd nts 


像 强 


Fig.7 Correlation between the power values (in arbitrary 
unit) received and the brightness temperature of the radio 
source after image intensity calibration at different 


frequencies on November 22, 2015 


图 8 展示 了 MUSER 观 测 到 的 宁 
和 其 他 结果 的 对 比 , 红色 菱形 符号 表示 MUSER 观 
测 到 的 值 , 红色 的 线 为 拟 合 观测 值得 到 的 曲线 , 黑 
色 的 线 表 示 的 是 Zirin 等 3 观测 的 结果 , 绿色 的 线 
则 是 Tan 等 t3 结 合 宁静 太阳 大 气 的 电子 密度 、 电 
子 温 度 模型 和 怀柔 频谱 仪 十 多 年 的 观测 数据 得 到 
的 结果 . MUSER 拟 合 曲线 介 于 两 个 其 他 结果 之 间 ， 
与 Tan 等 引得 到 的 亮 温 度 曲 线 具 有 非常 高 的 相关 
性 , 达到 0.96, 从 这 些 结果 也 反映 了 定 标 结果 的 准 
确 性 . 
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SPER 观 测 到 的 宁静 太阳 亮 温度 和 其 他 结果 的 对 比 


图 8 MU 


Fig.8 Comparison between the quiet solar brightness 


temperature observed by MUSER and other results 


4 结论 


在 射电 天 文 观测 中 , 射 MEL 度 是 天 体 
物理 研究 的 一 个 重要 物理 参数 ， 对 其 进行 定 标 是 射 
电 望 远 镜 数 据 处 理 的 重要 一 环 . 太阳 射电 观测 与 微 
弱 宇宙 天 体 的 射电 观测 不 同 , 由 于 太阳 射电 信和 号 较 
强 , 射电 日 像 仪 的 灵敏 度 要 求 并 不 是 太 高 . 太阳 射 
电 辐 射 的 动态 范围 很 大 , 且 其 物理 结构 比较 复杂 ， 
既 有 大 尺度 的 展 源 结构 , LAMI MEDI. HH 
等 相应 的 小 尺度 结构 . 因此 , 不 能 利用 传统 的 射电 
望远镜 定 标 方法 来 实现 对 射电 日 像 仪 的 图 像 强度 
定 标 . 
基于 太阳 射电 图 像 数据 中 含有 太阳 圆 盘 的 结 
构 信息 , 上 且 其 主要 反映 在 阵列 短 基 线 的 可 见 度数 据 
中 , 本 文 提出 了 一 种 适用 于 射电 日 像 仪 的 图 像 亮 温 
度 定 标的 方法 . 利用 短 基 线 的 幅度 信息 拟 合 第 一 类 
贝 塞 尔 函 数 , 可 以 得 到 图 像 中 宁静 太阳 圆 盘 的 射电 
半径 和 强度 . 再 利用 瑞 利 - 金 斯 定律 和 每 天 太阳 射 
电流 量 即 可 得 出 每 天 图 像 的 定 标 因子 G。, 从 而 实现 


6535 周志 超 等 : 明 安 图 射电 频谱 日 像 仪 图 像 的 强度 定 标 方法 3j 
对 射电 日 像 仪 的 图 像 定 标 . 6] Yan Y H, Zhang J, Wang W, et al. EM&P, 2009, 104: 97 

利用 文中 提出 的 方法 对 MUSER 观 测 的 太阳 图 7] Yan Y H, Chen Z J, Wang W, et al. FrASS, 2021, 8: 
像 进行 了 亮 温度 定 标 , 包括 宁静 太阳 和 有 太阳 射电 984043 


8] Thompson A R, Moran J M, Swenson G W. Interferom- 


j > A8 AS El Ip s y d Asi 42H H5 Se 

爆发 等 不 同 的 In. 结 末 显 示 每 天 得 到 的 定 标 因 si etry and Synthesis in Radio Astronomy. Third Edition. 
误差 较 小 ， 定 标 后 的 亮 温 度 与 单 天 线 接收 到 的 功率 New York: Springer Nature, 2017: 91-94 

值 变 化 趋势 相同 , 具有 很 强 的 相关 性 . 同时 , 得 到 不 9] ER, BBE, 刘 东 浩 , 等 . 电子 学 报 , 2019, 47: 4 


同 频率 IJ 37 B8 S EH 亮 温 度 与 其 他 文献 的 结果 也 比 10] Zhou Z C, Yan Y H, Chen L J, et al. RAA, 2022, 22: 
— 05019 

mU. iH etx E AEA nT E 

9 7 Hi 说 H 了 此 定 标 方法 的 有 效 PERUTI SE 11] Geng L H, Yan Y H, Chen Z J, et al. 2016 URSI Asia- 

Pacific Radio Science Conference, Seoul: IEEE, 2016: 


参考 文献 998-1999 


Ad bo á 4 " . 12] Zhang M H, Zhang Y, Yan Y H, et al. RAA, 2021, 21: 
1] BRE. 第 3 章 太 阳 射电 频谱 成 像 观测 // 史 建 魁 , rp zx, 刘振兴 . T T 


RE "m a Sia 284 
空间 物理 学 进展 (第 5 卷 ). 北京 : 科学 出 版 社 , 2015: 71-85 NS nr. A 
2] Aschwanden M. Physics of the Solar Corona: An Intro- I8] Airin H, Baumer B Ms Hutiord G J AD; T99779 


duction with Problems and Solutions. Second Edition. PH Fedotova A, Altyntsey A Rochanov dietal etal: SEP, 
Berlin: Springer, 2005: 639-640 2019, Quer 

3] Chen P F. LRSP 2011, 8: 1 15] Fleishman G D, Loukitcheva M A, Kopnina V Y, et al. 
ApJ, 2018, 867: 81 


FEE 


4] Lin R P. Particle Acceleration and Transport in the He- 


liosphere and Beyond. AIP Conference Proceedings, Ky- 16] Ranganath G S. Resonance, 2008, 13: 115 
oto, Japan, June 20-25, 2008 17] Penzias A A, Wilson R W. ApJ, 1965, 142: 419 
5] Chen X Y, Yan Y H, Tan B L, et al. ApJ, 2019, 878: 78 18] Tan C M, Yan Y H, Tan B L, et al. ApJ, 2015, 808: 61 


Intensity Calibration for the Mingantu Spectral Radioheliograph 
Images 


ZHOU Zhi-chao? CHEN Lin-jie! YAN Yi-hua'? WANG Wei! 


(1 State Key Laboratory of Space Weather, National Space Science Center, Chinese Academy of Sciences, Beijing 100190) 
(2 University of Chinese Academy of Sciences, Beijing 100049) 


AssrRAcCT The Mingantu Spectral Radioheliograph (MUSER), a new generation of solar dedicated 
radio imaging spectroscopic telescope, has realized high-time, high-angular, and high-frequency resolution 
imaging of the Sun in the 0.4-15 GHz ultra-broadband frequency range. The radio brightness temperature 
is an important parameter to describe the solar physical process. It plays a very important role in the study 
of different radio radiation mechanisms, solar magnetic field and the acceleration of non-thermal particles 
in the solar burst process. Therefore, the image of radioheliograph must be calibrated for brightness 
temperature. This paper introduces a method suitable for radioheliograph image intensity calibration. 
The solar radio image contains the structural information of the solar disk. T'he radio radius and intensity 
of the quiet solar disk in the image can be obtained by fitting the first kind of Bessel function with the 
visibilities of short baselines of the radioheliograph. Then, the Rayleigh-Jeans law and the daily solar radio 
flux can be used to calculate the calibration factor Ce of the daily image, so as to realize the calibration of 
the MUSER image intensity. Applying this method to the actual observation data of MUSER, including 
different situations such as the quiet sun and solar radio bursts, the error of the daily calibration factor 
Ge is within 1096 of its value, and the obtained bright temperature of the quiet sun has a high correlation 
with the result obtained by other people. These indicate the feasibility and effectiveness of this method. 


Key words Sun: radio radiation, instrumentation: radioheliograph, instrumentation: calibration 
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